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We study the linear magnetohydrodynamic (MHD) equations, both in the Newtonian and the 
general-relativistic limit, as regards a viscous magnetized fluid of finite conductivity and discuss in¬ 
stability criteria. In addition, we explore the excitation of cosmological perturbations in anisotropic 
spacetimes, in the presence of an ambient magnetic held. Acoustic, electromagnetic (e/m) and fast- 
magnetosonic modes, propagating normal to the magnetic field, can be excited, resulting in several 
implications of cosmological significance. 

PACS numbers: 


I. INTRODUCTION 


Although Plasma Physics and Cosmology are two well- 
established fields of Theoretical Physics, the formulation 
of magnetohydrodynamics in curved spacetime is a rel¬ 
atively new development [1]. In particular, in spite the 
fact that MHD processes in flat spacetime gained much 
attention [2] , only recently we were able to discuss exact 
spherically symmetric MHD solutions within the context 
of General Relativity (GR) [3], something that gave rise 
to efforts of exploring MHD processes in the vicinity of 
central engines [4], [5], [6]. On the other hand, magnetic 
fields are known to have a widespread presence in our 
Universe, being a common property of the intergalactic 
medium in galaxy clusters [7], while, reports on Faraday 
rotation imply magnetic fields of significant strength at 
high redshifts [8], [9]. These results lead the scientists 
to go even further and to look for gravitational instabili¬ 
ties in magnetized cosmological spacetimes, either in the 
Newtonian [10] or the GR limit [11], [12], [13], [14], 
[15], [16]. 

In the present article, we intend to re-construct the 
MHD equations, both in the Newtonian and the GR limit 
and to use them in order to study the finite-amplitude 
wave propagation in MHD media. The reason to do so, 
actually relies on the thermal history of the Universe. 

According to the Standard Model [17], [18], after nu¬ 
cleosynthesis the Universe goes on expanding and cooling 
but nothing of great interest takes place until t ~ 10^^ 
sec [z ~ 10^). At that time, the temperature drops to the 
point where electrons and nuclei can form stable atoms 
(recombination epoch). During the so-called radiation 
epoch, photons couple strongly with matter, the main 
constituent of which is in the form of plasma. Interactions 
between the various constituents of the Universal matter 
content along this period include radiation-plasma cou¬ 
pling which is described by plasma dynamics. On the 
other hand, the presence of plasma plays an important 
role in shaping the radiation spectrum, something that 
is fortified by the fact that there appear to exist cosmic 
magnetic fields, of the order 10“® G, which evidently 
have not had enough time to evolve during the gravitat¬ 


ing epoch of the Universe [7], [12]. Therefore, although 
it is not traditional to characterize the radiation epoch 
by the dominance of plasma interactions, it may well be 
also called the plasma epoch. It is the time-period during 
which the electromagnetic interaction dominates among 
the four fundamental forces. 


II. MHD EQUATIONS IN AN EXPANDING 
UNIVERSE 


A. The Newtonian MHD Theory 


During the plasma epoch, the MHD equations in the 
Newtonian theory of gravity are given by [19] [20] 


fDui 


dp dU 1 d{H,H^ - 
dxi ^ dxi 47r dxj 


1 Dp dui 

p dt dxi 


= -47rGp 


( 3 ) 


dH, _ djujHj - UjHj) dHj _ 
dt dxj ’ dxi 

In Eqs (1) - (4), Latin indices refer to the four¬ 
dimensional spacetime (in accordance Greek indices refer 
to the three-dimensional spatial section) , Sij is the Kro- 
necker symbol and we have used the notation 

D _ d d 
dt dt dxj 

On the other hand, p and p denote the energy-density 
and the pressure of a magnetized perfect fluid. Hi are 
the strength-components of the ambient magnetic field 
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and U is the gravitational potential. Introducing the 
scalars [21] 

(5) 


‘1(7 — (7ij(7ij , 


25" — (7ij(Jij , 


Ilu — uj(jj'ij 
loj — UJ'ijLO'ij 
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and 0 = is the expansion parameter, we obtain 


— (Jij -\- iJij -j- 


Combining Eqs (1) - (4) with Eq (8), we find 
1 D^p _ d ^Idp* 
dxi p dxi 


p df^ 

Hj dHi dp 
A-Kp^ dxj dxi 
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where, we have set p* = p + ^ to denote the general¬ 
ized pressure involved, due to the perfect fluid plus the 
magnetic field. Eq (9) governs the finite-amplitude wave 
propagation in Newtonian MHD theory. When = 0, it 
reduces to the corresponding law derived by Hunter [22] . 
On the other hand, using the notation 

d 


DH, .d ^d 

dt ~^dt^ d?^"~ dxj dxj 

together with Eq (8), we obtain 


DH, 


dt 


- = 


■UJ, 


ij 


( 10 ) 
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and multiplying Eq (11) by pHi, we end up an equation 
for the temporal evolution of the magnetic field’s energy 
density, namely 

D pH^ p AO pH^ 

where, p is the permeability [11]. 

Summarizing, the MHD dynamics in the Newtonian 
limit is determined by the following system 

,Dui dp dU 1 d{HiHj — 

P^-;r) = -TGr + PTGr + 


dt ' dxi dxi At: 
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B. The Relativistic MHD Theory 

To derive the general relativistic MHD equations, we 
start up with the Einstein field equations [19] 

1 SttG , , 

Hah T^dabH — idTah ; ^ — 7 (17) 

2 ed 

where, Hat is the Ricci tensor, TZ = gab'H°‘^ is the 
corresponding scalar curvature and Tab is the energy- 
momentum tensor, representing the average state of 
the matter-energy content. Eq (17), together with the 
Bianchi identities 

(T^ab _ = 0, (18) 

result in the conservation law T°‘\ = 0, where, the semi¬ 
colon denotes covariant derivative. In what follows, we 
perform our calculations in the system of geometrical 
units, where G = c = 1, admitting that the fluid’s four- 
velocity satisfies the condition Uau“ = —1; i.e. we place 
ourselves in a coordinate system comoving with the fluid. 

Again, we consider a magnetized perfect fluid source, 
determined by the following energy-momentum tensor 

r“^ = (p + :|!)MV + (p+^)/i“^-iJ“ff^ (19) 

where. 


f^ab ^ gab ^a^b 

is the projection tensor and 

e = p -I- pH 

is the total energy density, due to the mass-energy con¬ 
tent and the internal motions (pH). By analogy to 
the corresponding evaluation in the Newtonian limit, 
we introduce the expansion 0 = rtfjj, the shear dab = 

K^b'^ic-,d) - \Phah and the twist ujab = K^b'^lc-.d] (round 
brackets denote symmetrization while square brackets 
antisymmetrization), to write the covariant derivative of 
the fluid’s four-velocity in the form 

Ua-,b = (Jab + bJab + ^dhab “ UaUb (20) 

where, we have set = u“ u'”. 

Accordingly, the MHD dynamics in the GR limit is 
determined by the following system [11] 

tt2 

xu‘^ + -k a;0M“ + {p+ —= (77“i7^);b (21) 

iP - ^labU'^U^ = h‘^\p + ^),ab + 2^(ff'0) - 

— {H‘^H^)-ab + 2^(“^-h cr2 — — iHiia) T 

|(p + 3p + H^) + 2Ua{H-H%b + (i7');aU“ 


(16) 
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where, we have set x = p + p + H^. Eqs (21) - (23) 
may be applied to the investigation of perturbation ef¬ 
fects in curved spacetime and, hence, to the discussion of 
linearized stability criteria. In order to obtain the per¬ 
turbed expression of the above equations, we follow the 
well-established method imposed by Hawking [23]. 


III. PART I: MHD PHENOMENA ALONG THE 
PLASMA EPOCH 

To discuss MHD phenomena in the early stages of 
the Universe, we use a solution to the Einstein-Maxwell 
equations, first developed by Thorne [24] and later dis¬ 
cussed by Jacobs [25], which contains an ambient mag¬ 
netic field. The model is homogeneous and has two equiv¬ 
alent ’’transverse” directions and one inequivalent ’’lon¬ 
gitudinal” direction at each point of the spacetime. The 
magnetic field is frozen into the matter content and is di¬ 
rected along the longitudinal direction (z-direction). We 
consider that the matter-content filling this model is in 
the form of the perfect fluid (19), which obeys the equa¬ 
tion of state p = 7p, with 1/3 < 7 < 1. In this case, the 
line-element reads [24] 

ds^ =—dt^+t{dx'^ + dy'^)+ t'^dz'^ (24) 


Upon consideration of the perturbed MHD equa¬ 
tions [23], after straightforward but tedious calculations, 
we obtain the equation governing the density fluctua¬ 
tions [11] 


= ci5u° + + FiSH^ + F 26 H^^f, 


■ ~ [^3 + Cs^4]i5p 
+ FaJi/ffca (26) 


where, = Sp/Sp is the speed of sound and Ai,Ci,Fi 
{i = 1, 2, 3), are all functions of t and 7, the explicit form 
of which is given in [11]. Inserting Eqs (25) into Eq (26), 
we obtain a secular equation of the form 

— inA — Si — ikS 2 = 0 (27) 


where. A, Si and S'2 are complicated expressions of the 
background energy density and the pressure, as well as 
of the shear and the expansion of the Universe, together 
with the magnetic field’s strength (see Appendix A). Eq 
(27) results in 

ni = r cosw/2 -I- i[— -I- r sinw/2] (28) 

and 

712 =—Tcosta/2 -I- i[— — r sinw/2] (29) 

where, we have set 


where, m = 2p^p^. With the aid of Eq (24), we may 
apply the GR limit of MHD dynamics to discuss sta¬ 
bility of Bianchi-Type I cosmological models. In this 
case, we depart from the traditional approach, adopting 
an anisotropic cosmological model, in which the inherent 
magnetic field is responsible for the initial anisotropy of 
the Universe. 


A. Excitation of low-frequency plasma waves 

To deal with small-amplitude waves, we first assume a 
background situation, representing a uniform plasma in 
curved spacetime, in the presence of an ambient magnetic 
field. In our case, this situation is described by the so¬ 
lution (24), together with the energy-momentum tensor 
(19). Accordingly, we introduce first-order perturbations 
to the MHD equations (21) - (23) and taking into account 
the so-called Cowling approach [26], which admits that 
dpab = 0, we neglect all terms higher or equal than the 
second order. Eventually, we search for solutions to the 
linearized MHD equations in which all the perturbation 
quantities are written in the form of a plane-wave, with 
angular frequency n and a wave-vector k = (0, 0, k^) par¬ 
allel to the magnetic field, H = (0,0,77^). Accordingly, 

(25) 


and 



AF 


k^Sl 




4Si - A^ 
4t 


(30) 


(31) 


According to Eqs (28) and (29), the medium remains 
stable when the frequency equals to rii. In this case, 
the magnetic field’s perturbation, turns out to be 

an Alfven wave. On the other hand, for n = n 2 , the 
medium becomes unstable and SF[^ grows in time. We 
identify this instability with the well known Parker’s in¬ 
stability [27]. 

There are also other potential instabilities that can be 
traced within the context of general relativistic MHD. In 
particular, a spectrum of magnetized sound waves may be 
excited and form large-scale damped oscillations in the 
expanding Universe (e.g. see [16]). The characteristic 
frequency of the excited waves is slightly shifted away 
from the sound frequency and the shift depends on the 
strength of the primordial magnetic field. This magnetic- 
field-dependent shift, may have an effect on the acoustic 
peaks of the CMRB. 

One of the main scopes of this study, is to exam¬ 
ine the spectrum of the unstable low frequency plasma 
waves in the anisotropic cosmological model (24). Under 
the assumption that the magnetic field lays along the 


(5e, 5p, SF[°', 5u°' ~ expi(7it — kz). 
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z-direction and the perturbations are plane-waves prop¬ 
agating also in the i-direction, the perturbed general rel¬ 
ativistic MHD equations, result in the dispersion relation 


^ ^2)[__02 + + -(p + 3p) 

2eHl 46»2ij2 2ij2 1 


— iri[40 -I- 


2^,0 


40H^ 


ip + p) Hp + p) 

= {TZi -I- inTi2)^^— 
de 


(32) 


where, we have set 

SH^ 




Se 47 rl ?(2 — 7 ) 


[1 — Ai — inA2], 


(33) 


and 72.1, 72-2, Ai, A 2 and V, are functions of the magnetic 
strength H, the shear a and the expansion 9 of the back¬ 
ground solution (see Appendix B). At the early stages of 
the evolution {t 0), Eq (32) reads 


[1 - -f -k -k Ti ( 7 ,7) -f inT2{t, 7) = 0 

(34) 

where, ua is the Alfven velocity and 11 ( 7 , t), T 2 {'y, t) are 
complicated expressions of 7 and 7, the exact form of 
which is given in Appendix B. We re-write Eq (34) in 
the form Dr + iDi = 0 and, assuming that n = Ur + iui, 
where, the real part of the excited frequency is much 
larger than the corresponding imaginary part (n^ >> 
rii), we obtain: 


• The real part of the excited frequency from the 
equation Dr = 0 


Ua , ,22 k ^ u \ "^1(7) 


(35) 


• The imaginary part of the excited frequency from 
the relation [28] 


77j (tZj. , k'j 

dDr(nnk) \ 

an In=nr 


which yields 


2(1-u^(7)/(2-7)) 

= fm] ( _^^ 


(36) 


(37) 


Eor t ^ 0, Ui becomes negative, decaying rapidly 
at late times. Therefore, we verify that the role of 
magnetic field is to make the expanding Universe 
even more stable against the expected Jeans-type 
instabilities. 


On the other hand, the real part of the frequency is 
shifted away from the frequency of sound 
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2 
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-I- 


2-7 
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[it) 


2-7 


II 

*2 


(38) 


In fact, as 7 ^ 1/3, the magnetic field vanishes 
and the anisotropic model (24) approaches the weakly- 
magnetized Friedmann - Robertson - Walker (FRW) 
model, which is used extensively in the literature (e.g. 
see [31] and references therein). In addition, the char¬ 
acteristic frequency approaches the value Ur = kcg -\- 
Anr{uA). The presence of the Alfven velocity in Eq (38), 
may be responsible for the distortion of acoustic peaks, 
as it was pointed out first by [12]. Therefore, a spectrum 
of low frequencies will be excited with frequency Ur . 


B. Evolution of the density fluctuations 

Finally, we examine the temporal evolution of the den¬ 
sity fluctuations. A linear density perturbation has the 
form 


Sp 


Poe 


\Bil,t)\ 


^i{nrt — kz) 


(39) 


where. 


- 8 ( 7 , 7 ) 


g(7) 

2-7 


The function j 77 ( 7 , 7 )] decays as time grows or, equiva¬ 
lently, as 7 decreases. In particular, as 7 —> 0, which cor¬ 
responds to 7 —> 1, we have ua ^ 1 and | 7 ?( 7 , 7 )j —> 00 . 
On the other hand, for 7 —> 00 , which results in 7 —> 1/3, 
we obtain ua —> 0 and l7?(l/3, 7)j —> (3/2). These results 
indicate that, at early times (t 0 ), the perturbations’ 
amplitude becomes large and the excited waves will form 
a spectrum of damped oseillations within the magnetized 
cosmological model. We may search for the source of this 
damping by analyzing the functional form of T 2 {'y,t) in 
Eq (34). In the absence of magnetic fields (7 = 1/3), 
we obtain T 2 = —40 and the expanding plasma stabilizes 
the ion-acoustic waves. 

On the other hand, if the cosmological model is static 
{A = const, and 0 ~ ^ =0), stabilization of the mag- 
netosonic waves can be caused by the magnetic pressure. 
In a magnetized expanding Universe [Eq (24)] both pro¬ 
cesses will be combined to damp the large-amplitude os¬ 
cillations. Assuming, for the sake of agreement, that at 
a some time 7 = const, the Universe becomes static, a 
Jeans-like instability sets in (see also [29]). Therefore, 
we realize that the damping is caused by the free energy, 
available both in the plasma and the magnetic field, dur¬ 
ing the early stages of the Universal expansion. 

Summarizing, we have analyzed the small-amplitude 
wave propagation in the early Universe, using an 
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anisotropic cosmological model with an inherent mag¬ 
netic field towards the z—direction, which is frozen into 
a perfect fluid with equation of state p = (where, 
1/3 < 7 < !)■ For 7 —> 1/3, the anisotropic model 
used, approaches the weakly magnetized FRW model, in 
which the magnetic field is a small linear perturbation 
(e.g. see [14], [30]). 


V • B = 0 (44) 

where, a(t) is the radius of the Newtonian Universe and 
(j), stands for the gravitational potential. 


IV. PART II: MHD PHENOMENA OF 
COSMOLOGICAL SIGNIFICANCE 


Newtonian Cosmology can be a quite good approxima¬ 
tion for gravity, at scales much smaller than the Hubble 
length. For this reason, it is often used to describe our 
Universal neighborhood (0.1 < z < 1.5), within the con¬ 
text of Observational Cosmlogy; in particular, as regards 
the local observations of cosmological significance [18]. 
Within the context of Newtonian Cosmology, we study 
the linear MHD behavior of a viscous magnetized fluid 
with finite conductivity and generalize the Jeans instabil¬ 
ity criterion. In particular, we look into the linear evolu¬ 
tion of small inhomogeneities in the cosmic medium and 
examine how the magnetic field and the fluid’s viscosity 
may affect the characteristic scales of the gravitational 
instabilities. Accordingly, we discuss the electrodynamic 
properties of the collapsing fluid, the resulting amplifica¬ 
tion of the magnetic field and the formation of unstable 
current-sheets. 

A central point to our discussion is the concept of 
anomalous resistivity, which is triggered by electrostatic 
instabilities in the plasma and may reduce substan¬ 
tially the corresponding electrical conductivity. We argue 
that these changes in the resistivity of the protogalactic 
medium might have lead to the formation of strong elec¬ 
tric fields during galactic collapse. These fields can ac¬ 
celerate the abundant free electrons and ions, to produce 
Ultra High Energy Cosmic Rays (UHECRs) during the 
formation of protogalactic structures. 

In Newtonian Cosmology, one deals with the following 
set of MHD equations [19] 


dp 

m 


„ a 1 ^ 

—3—p-V • (pu) 

a a 


(40) 


du 

'm 


-h 


o- ^ c, ^ 1^, 

—u -(u-V)-u-Vp -I —Vd 

a a ap a 

—(V X B) X B -\- 
47rap a-^p 


(41) 


= -47rGaV (42) 


dB 

dt 



(43) 


A. Magnetically induced anisotropies in the 
gravitational collapse 


Once again, we perturb Eqs (40) - (44) around a ze¬ 
roth order background solution, representing a uniform 
plasma in the presence of an ambient magnetic field 

/O = po + Pi , B = Bo -\- Bi , (j) = (po -\- (pi , u^O 


We assume plane-wave perturbations and take the time- 
derivative of Eq (41). The real part yields 


= - H 


vk^ 
a^po 

„ k^, o 9 vH , 
SttGpo-9(^8 + '^A 4 -) 

Po 


(45) 


Furthermore, we assume that the perturbations propa¬ 
gate (a) parallel to the magnetic field and (b) perpendic¬ 
ular to the magnetic field. Accordingly, we end up with 
two critical wavelengths, associated to the corresponding 
Jeans scales, namely 


Aj. ~ 


lcj-\-u\-\- vH/po 


SttGpo 


All ~ 


I -\- vH/po 
SttGpo 


(46) 


where, H = ^ is the Hubble parameter, v is the viscosity 
of the plasma fluid and p is the corresponding resistivity. 
From Eq (46), it is evident that the magnetic field induces 
a degree of anisotropy in the gravitational collapse and, 
therefore, one expects gradual formation of turbulent mo¬ 
tions within the magnetized medium. These motions, in 
turn, contribute additional terms to the viscosity of the 
fluid. In the post-recombination epoch, the viscosity due 
to turbulent motions reads 


^turb ^ PO'^llmix (47) 

where, Imix the turbulent mixing length. Assuming that 

• ui reaches values close to ua 

• Imix A^ and 

• Tree is known 

we end up with the condition 


cl < 


vH 

PO 


<< 


(48) 
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At the present epoch, we adopt the typical values [18]: 
For the magnetic field, B < 10“^ G; for the energy- 
density, po > 10“^® gr • cm“^ and for the Hubble param¬ 
eter, H = 100 ft. km • sec“^ • Mpc“^, where 0.4 < ft < 1, 
we find 


All 


' cIpo 

^turhH 


> 1 


(49) 


with Ax being of the order of a (comoving) Mpc. Follow¬ 
ing the standard structure formation scenarios, the initial 
collapse of a large structure will be followed by successive 
fragmentation into smaller scale formations with charac¬ 
teristic lengths A ^ Ax. Moreover, as we will outline 
next, the magnetic field which is trapped into the gravi¬ 
tating medium, will be increased even further, due to the 
anisotropic collapse. In the case of an almost spherically 
symmetric collapse, linear inhomogeneities in the mag¬ 
netic field’s energy density amplify in tune with those in 
the density of the matter, so that 


SB^ oc Sp 


where, 6B^ = B^/Bq and Sp = pi/po Therefore, even in 
the spherically symmetric case, the formation of matter 
condensations in the post-recombination Universe signals 
the amplification of any magnetic field that happens to 
be present at that time [32]. 

We have already seen that the anisotropic magnetic 
field induces some degree of anisotropy in the gravita¬ 
tional collapse of the protogalactic structure. On the 
other hand, this magnetically induced anisotropy will re¬ 
act back to affect the evolution of the field itself. The evo¬ 
lution of the magnetic field during the nonlinear regime of 
a generic, non-spherical protogalactic collapse, has been 
considered by a number of authors (e.g. see [31] and 
references there in). The approach is both analytic and 
numerical and agree that shear effects increase the final 
strength of the magnetic field, while confining it onto the 
protogalactic plane. Compared to the magnetic strength 
of the spherical collapse scenario, this anisotropic in¬ 
crease is larger by, at least, one order of magnitude and, 
therefore, protogalactic structures can be endowed with 
magnetic fields considerably stronger than those previ¬ 
ously anticipated. 

So far, we have seen how the magnetic presence mod¬ 
ifies the way the gravitational collapse may proceed, by 
changing the overall stability of the magnetized fluid. 
This, in turn, affects the evolution of the magnetic held 
itself and may trigger a chain of nonlinear effects on cer¬ 
tain scales. 

At next, we will see that, selective amplification of cer¬ 
tain perturbation modes, plays an important role during 
the nonlinear stages of protogalactic collapse, helping in¬ 
stability to reach its saturation point. 


B. Amplification of the induced electric currents 


In the presence of magnetic fields there is an electric 
current induced by B, namely 

J=-^VxB = -^aB (50) 

Att Att 

where, V x i? = aB and a measures the magnetic tor¬ 
sion (see [33] for details). One expects that, initially, a 
is small. However, the subsequent fragmentation of the 
protogalactic cloud to scales with A ^ Ax will increase 
V X H and strengthen the induced current. In particular, 
for B ~ 10“^ G, we obtain J ~ 10®a which may reach 
appreciable strengths for reasonable values of a. 

When the induced current exceeds a critical value, 
namely Jc ~ p{elmp)cs, where, Cg ~ I0'^y/T{°K) 
cm/sec, rup is the ion mass and e is the electron charge, 
the excitation of low frequency electrostatic turbulence 
will increase the resistivity of the medium by several or¬ 
ders of magnitude [34]. In fact, even small values of a 
may lead to J > Jc, thus making the plasma electrostat¬ 
ically unstable. This effect, which is known as anoma¬ 
lous resistivity, can be explained by the development of 
current-driven electrostatic instabilities in plasma, that 
could lead to the excitation of a variety of waves and 
oscillations. Absorption of these waves by the ions, re¬ 
sults in a transfer of momentum from the electrons to the 
ions, along with the usual momentum loss from the for¬ 
mer species to the latter. The average momentum loss by 
the electron per unit time can be written as an effective 
collision term, in the form 


nrUeVeffUe = -Ffr 


(51) 


where, Ffr is the average friction force and n = p/rrip is 
the ambient number density of the plasma particles. The 
friction force is proportional to both the linear growth 
rate of the electrostatic waves (yj,) and the energy of 
the excited waves (H4). The effective collision frequency 
is estimated to be Vcff = oJefWsat/kBT). Then, the 
anomalous resistivity will be 


'Han 




ksT J V47ry We 


(52) 


where, Wg = 5.6 x sec“^ is the plasma frequency, 

fts is the Boltzman constant and {Wsat/ksT) ~ 10“® — 
10“^ is the saturated level of the electrostatic waves [34]. 
For certain types of current-driven waves, the anomalous 
resistivity is several orders of magnitude above the clas¬ 
sical one, as confirmed in numerous laboratory experi¬ 
ments [35]. 

This sudden switch to high electrical resistivity leads to 
the formation of strong electric currents and therefore to 
a fast magnetic dissipation and intense plasma heating. 
The electric fields induced by the gravitational collapse 
will be 

E ~ —B -t- rjanJc ~ TlanJc (53) 

c 
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provided that ^ ^ 1. According to this scenario, the 
gravitational collapse of a magnetized dust cloud ampli¬ 
fies the magnetic field and (indirectly) generates strong 
electric currents, which form large-scale current-sheets 
localized on the protogalactic plane. 

The anisotropy of the collapse enhances the local cur¬ 
rents even further and, eventually, drives the resistivity 
towards anomalously high values. The inevitable result is 
strong electric fields accelerating the abundant free elec¬ 
trons. In particular, the energy gained by an electron 
travelling a length, A ~ 10“^A_l Aj_ is Wkin ~ eAA 
and the relativistic factor, 7 = [1 — (v/c)'^]~^ is given by 


Wkin ^ eAA_L ^ lQe^r]anVTB 
meC^ lO^meC^ mempC^G^/^ 



in CGS units. Taking into account that 
Ax ~ B/niJipVG, 

Cs < ^H/p < u\ 
and, through Eq (52), 

Wsat/ksT 

we find that the typical energy gain by a free electron 
may reach extremely high values within short timescales 
{tacc ~ A/c ^ Ax/c ~ 10® yrs), even for relatively weak 
magnetic fields. It is straightforward to extend this pro¬ 
cess to proton acceleration and show that protogalactic 
collapse may be the source of UHECRs. 

We also anticipate a few particles travelling across sev¬ 
eral fragments adding up to a scale comparable to the 
Jeans length. Many particles will drift in and out the 
primordial current-sheet, as well as the associated strong 
electric fields. Through synchrotron emission, particles 
gain and lose energy continuously. Provided that frag¬ 
mentation has already taken place, these particles may 
diffuse along the different current-sheets, forming the 
observed power-law distribution. The role of unstable 
current-sheets on the acceleration of cosmic rays along 
the giant radio-galaxies, has already been pointed out in 
the literature [36]. 

Summarizing, the suggestion made here is that, in the 
context of Newtonian Cosmology, production of UHE¬ 
CRs may have started almost simultaneously with the 
formation of galaxies, through the electrodynamic char¬ 
acteristics of the gravitational instability. This process 
may be continued until today, since the previously de¬ 
scribed instability is active on all cosmic scales. 


V. CONCLUSIONS 

Part I: 

1 The Jeans instability, established in an Einstein 
Universe with a very particular cosmological con¬ 
stant [29], is absent in the presence of magnetic 
fields, which act as stabilizing forces. 

2 An anisotropic magnetized Universe re-enforces the 
stability, initially found in the FRW model, against 
Jeans-like instabilities. 

3 A spectrum of low frequency and large amplitude 
damped oscillations may appear in the early Uni¬ 
verse, with a characteristic frequency = kc^ + 
An, where. An depends on the strength of the pri¬ 
mordial magnetic field. 

Because of the above results, we come to the following 
conclusions: 

• The magnetic-field-dependent shift (An), may 
cause measurable distortion on the acoustic peaks 
of the CMRB (see also [12], [15], [30]). 

• The amplitude of the excited waves is substantially 
large in the early stages of the Universal evolution 
and gradually decays. 

• If the Universe has passed through a strongly mag¬ 
netized anisotropic phase, the excited waves may 
be responsible for the formation of large-scale fluc¬ 
tuations, as it has been shown in 2-D numerical 
simulations [37]. 

Part II: 

1 We have found that, a magnetic field frozen onto 
a fluid of finite viscosity, may alter the standard 
picture of the gravitational instability. 

2 We have discussed how the preferential, anisotropic 
magnetic amplification may increase the induced 
electric currents on a plane perpendicular to the 
main axis of the collapse. 

3 These gravitationally induced current-sheets will, 
in turn, trigger electrostatic instabilities, which 
may lead to anomalously high values of the resis¬ 
tivity and, consequently, to strong electric fields. 

4 These electric fields can be strong enough, in order 
to accelerate free electrons up to ultra high ener¬ 
gies, producing UHECRs. 

To the best of our knowledge, this is the first time 
that a direct connection between gravitational instabil¬ 
ity and cosmic-ray acceleration has been suggested and 
discussed. 
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Appendix A 


The expressions A, and S 2 , appearing in Eq (27), 
are written in the form: 


A = 


P + P 3 


I FT \ + 2p + , 


7^l = [773.00 - 4(0773.0 + 0oi?3) - 2i73r"3] 

+ 773[—-77^ +-0^ + 2cr^ + 2(p + 2p)] 

+ /"[-2r°3_o-2r03r^o + 2r°3ry773 
- in^ + e)H3 

7^2 = 2773,0 - 40773 - r°37735'" 

Ai = ^[fc^(773)2 + (27-l)0[2xrio + (p + 77V2).o] 

- k^dg^^e+p)] 


Si — ~^(1 + Cs)(l + 3cs) 

,,, , H^-(p+P + 3H^)H 

+ ip + p+H^) J 

- clH^ - ^ [P + 3p + ^ + 4(7^ + ^772 - ^77^77,;o 


+ 


20 

P + P 
1 


(77,77.' 


.0 ^ ^0H^ + y077,77-) 


+ 2^^^3zz-,oHzi‘2:p + 2p + H^){p + p + H^) ^](1 + Cs) 

- {77foo - 6077.^) - 77^ [3772 + 6 a ^ - 60^ + 

O 

1 ,, zz /i\ 4“ 2p + 77^, 

+ nff 9zz,o(nff 9zz,0 ~ ~ . TJ 2 ] 

2 2 p + P + 77^ 


8 




;0^z;0 -I- HzH^qq + OHzH^q 


+ 2H^[2a^ + ^ + ^{p + p + Hd] 

' -(77,77."o + Hpi/2 ^ ) 


P + P 


277,77fo 


yrTTfp 


A 2 = (27 - l)a ;0 

V = -n 2 x + fc 2 ( 773)2 - ( 2 - 7)0 


22;r3Q + {p+ ~).o 


77 


2a:rgo + p + — - (2 - 7)x0 


Finally, the functions Ti and T 2 , appearing in Eq (34), 
have the following functional form: 


ri(7,7) = 


TM. 


fi(7) = -(97®-9937®+ 48427^- 68187^+ 28897^ 

- 4617 + 788)7(6(7- 2)(7- 3)2(1+7)2) 


and 

T2dA) = -40 + 


40 / 


772 




3 \47r(p + p) 
2773.0773 


40 


772 


47 r (2 — 7 ) 


772 


+ p""5...o [57720 _ ^^pjz _ ^2^z ^ , 

P + P + 


47 r(p + p) 47 rp (2 — 7 ) 

(1-7) ( 


jj2 

~);0 ' (7 + l)t V 47 rp (2 — 7 ) 


p + p + H'^ P + P 


S 2 = 77,77.^0 + - 57720 _ 


772 


IZZZ.Q 


P + P + H'^ 


(P- 


772 


where, 0 = —2/[(l + 7 )^]. The function T 2 depends on 
7 and t, since the magnetic Held, the energy density and 
the pressure are all functions of 7 and t. In other words, 
T 2 ( 7 ,t) has the form 


+ g^^H 


2 2 p + 2p + 772 

p + p + 772 


12 ( 7 , 7 ) = 


C(7) 


where. 


Appendix B 


C(7) = 


457^ - 367^ - 12672 + 2127 - 127 

3(7-3)(7-2)(7+1)2 


The functions 77i, 7 ^ 2 , Ai, A 2 and V, appearing in Eq We can easily show that ri( 7 , t) is positive and T2(7, *) 


(32) are: 


is negative, for all the allowed values ofl/ 3 < 7 <l. 
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